Optical conductivity of charge carriers interacting with a two-level systems reservoir 
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I. INTRODUCTION 

The possibiHty of manipulating physical systems at 
length scales where quantum effects become important 
has attracted much attention over the last years. This 
scenario has stimulated considerable effort in the study 
of quantum open systems, firstly because of its relevance 
to the phenomenon of quantum coherence and secondly 
because it might not be entirely obvious the way in which 
we can obtain information about the quantum dynam- 
ics of the system of interest coupled to its environment. 
One of the strategies broadly used in solving this prob- 
lem consists in the replacement of the environment by an 
approximate model. This must be done in such a way 
that after tracing the environment coordinates out, the 
problem can be formulated only in terms of the variables 
of the system of interest ji It is remarkable that most of 
the environments can be represented either in terms of 
a bath of oscillators2»4, when the physics is dominated 
by the delocalized modes, or by an spin bath^, when the 
localized modes play the major role. In this paper we 
will be interested in another type of thermal bath that 
can be thought of as the projection of the usual oscilla- 
tor modes onto their two lowest lying levels. At very low 
temperatures, these truncated two-level systems (TLSs) 
have the same properties as the usual harmonic oscilla- 
tors, that is, the two baths exhibit the same quantum 
limit. However, as we will show, they strongly differ in 
the classical regime, at high temperatures. 

It is now well known that when the oscillator model 
with linear spectral density is used to mimic a ther- 
mal bath interacting with a quantum particle, the wave 
packet associated with the latter undergoes a damped 
motion, exactly as in the classical problemi^i In this situ- 



ation, and within the long time approximation, the aver- 
age over the environment variables results in an equation 
of motion for the particle without memory effects. There- 
fore, it directly follows that the transport properties of 
the quantum particle can be simply described in terms of 
the damping and diffusion coefRcients. As a consequence, 
the optical conductivity o-(ijj) of a single particle coupled 
to an oscillator bath has, in the so-called ohmic case, only 
an incoherent part, which has a Drude-like form. In this 
case the Lorentzian width is determined by a temperature 
independent damping constant. 

The phenomenological approach of representing the en- 
vironment by an oscillator bath was successfully used 
in the study of dissipative effects in macroscopic quan- 
tum coherence^ (the spin-boson model) and macroscopic 
quantum tunneling^. Moreover, there are particular sit- 
uations, see 1^ for instance, in which the oscillator model 
can be derived from microscopic theories just following 
the prescriptions of Feynman and Vernoni. That is also 
the case of solitons, whose transport properties can be 
investigated using the collective coordinate quantization 
schemeiAflr. In those cases the effective equation of mo- 
tion for the center of mass of the soliton leads, in the 
long time regime, to a temperature dependent damping 
constant. The form of the damping constant is such that 
the optical conductivity of a system of non-interacting 
solitons has again a Drude-like form, and in the low tem- 
perature limit correctly reproduces the finite free par- 
ticle Drude weight at zero frequency. This behavior is 
completely general for solitons and therefore indepen- 
dent of the non-linear field theory that supports these 
localized solutions. From |3| and we can conclude 
that, within the long time regime, the optical proper- 
ties of quasi-particles coupled to an oscillator bath have 
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always the trivial Drude form. Therefore, if results at 
variance with the latter are obtained in a{uj) measure- 
ments, they can not be attributed solely to the above- 
mentioned particle-reservoir interaction. In those cases 
more complicated models are requiredifi or alternative 
thermal bath descriptions should be employedi"^ 

In this paper, our main goal is to present a different 
kind of thermal reservoir which, in addition to a dissi- 
pative dynamics for the charge carriers, induces a non- 
trivial optical conductivity. Our starting point will be the 
problem of a single particle subject to a complex potential 
which can be represented as a distribution of local two- 
level systems. Using the well known Feynman- Vernon 
formalisroi we are able to traced out the environment 
modes and obtain the reduced density operator of the 
particles. This operator will be expressed in terms of an 
effective action containing all the information about the 
coupling of the carriers to the reservoir and provides us 
with an equation of motion for the particles from which 
the transport properties can be computed. In this way it 
will be shown that the optical conductivity of this system 
has a temperature dependent non-Drude behavior with 
a very rich structure which contrasts with the oscillator 
bath transport properties. 

This paper will be divided as follows. In Sec. II we 
present a model describing a single particle interacting 
with a set of TLSs. In Sec. Ill an effective equation of 
motion describing the particle dynamics will be derived 
and in Sec. IV the optical conductivity of the system 
is calculated considering different characteristics of the 
thermal reservoir. Finally, in Sec. V, we discuss our re- 
sults and present the conclusions. 



II. THE MODEL 

As it was already mentioned, we will be interested in 
the transport properties of a single particle coupled to a 
generic TLSs reservoir. This problem can be completely 
described by the hamiltonian 

H = Ho + Hr + H^, (1) 

where Hq stands for a particle placed in an external elec- 
tric field, and it is given by 

Ho^ — + exE. (2) 

The distribution of TLSs playing the role of thermal 
reservoir in (QJ is denoted by Hr and may be expressed 
in the A;— space as 

Hr = ^—a^k, (3) 

fc=i 

where Uzk is the z-Pauli matrix. Finally, the interaction 
between the particle of interest and the thermal bath Hi 



v(x) 




/-I / 



Figure 1: Complex potential with multiple local quartic struc- 
ture. The effective TLS at site / is described in terms of the 
detuning e between the wells ground states and a typical ma- 
trix element for the tunneling processes A. 



will be taken in such a way to induce transitions between 
the states of each TLS. A suitable choice is 

N 

Hi = -x^Jk(Jxk, (4) 

k=l 

where Jk and axk are, respectively, the coupling constant 
and the x-Pauli matrix. 

Although the problem defined by Eqs. (QJ-Q was not 
derived from a microscopic description of a concrete phys- 
ical system, it is still very useful because there are many 
systems which, under certain circumstances, behave as a 
truncated TLS. Indeed, in the process of modeling defects 
in crystalline solids or amorphous materials, one has to 
deal with a distribution of locally quadratic-plus-quartic 
potentials, as it is shown in Fig. In fact, if each local 
double-well in the distribution has fairly separated min- 
ima, and the typical energy scale huc (obtained by the 
harmonic approximation about each minimum and as- 
suming they do not differ much) is such that fkuc » kT 
(see Refs. ^ and ,11] for details) all the locally quartic 
potentials at positions I can be effectively described as 
two-level systems and 

Hi = -^{Aiaxi - siazi) , (5) 

where A; is a typical matrix element for the tunneling 
process between the two minima and ei is the "detuning" 
between the ground states in the two wells. This fact 
suggests a possible phenomenological description for a 
particle moving in a complex potential as shown in Fig. 1, 
where the source of dissipation is the induced transition 
between the states of the TLSs ensemble and may have 
the form 1^. However, in the course of this procedure one 
should keep in mind that although the rate Si/kT may 
have any value, the TLSs representation is valid only if 
hwc » kT ~ (?iA/,e/). Therefore the high temperature 
limit should be understood as hwc >> kT > Ml where 
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n is the cutofl frequency of the system, which can be 
assumed to be of the order of max(?iA;, e;). 

Knowing the details of the model describing the TLSs 
array we can go further and study the dynamics of a 
particle coupled to this environment. 



In the expression above 

So[z] = I 
Jo 



^z\t') + ez{t')E{t') 



dt' 



III. EFFECTIVE PARTICLE DYNAMICS 
A. The Influence Functional 

This section is devoted to investigating the effective 
dynamics of a particle interacting with the TLS thermal 
reservoir. We will use the Feynman- Vernon functional 
integral approach, which begins with the calculation of 
the reduced density operator of the particle of interest, 
namely 



p{x,y,t) = Tr[{x\e-''- p{0)e'-\y)] 



(6) 



where we are using the coordinate representation for the 
particle states, Tr denotes the trace over the bath modes 
and H is the Hamiltonian of the total system given by 

The density operator of the whole system at time t — 
will be assumed to be separable, p(0) — pp{0)pr{0), where 
p and r denote the particle and the reservoir, respectively. 
The reduced density operator ^ can be written as 

P{x,y,t) ^ / dx' / dy'pp{x',y',0)j{x,y,t;x'y',0), 



where the superpropagator J has the form 

J= f Vx{t') r Vy{t')ei^^'>^''^-'^°^y^^T[x,y]. (7) 

J x' J y' 

I 



corresponds to the action of a free particle in the presence 
of an electric field and J- denotes the infiuence functional, 
which is given by 



p,(0) [Tcxpj: I HMn)dt' 



(8) 



Texp~-J H,{x{t'))dt' 



where T indicates time-ordered product and 

H,{z)^e^"'''H,(z{t'))e'*"''\ 



The central quantity in this method is the influence 
functional defined in Eq. Our next step is to de- 

rive an explicit expression for it. Although this quantity 
has already been evaluated through different approaches 
for interactions of the form I^JJ^, here we will sketch its 
derivation for the sake of completeness. 

In order to proceed, we will assume that the interaction 
strength is weak enough, such that we may expand Eq. 
in powers of Jk and retain only terms up to second 
order. We then obtain 



T[x,y] = ^-^l dt' dt" [(HMt'))HMt"))) + (HMn)HMt'))) - {H,{y{t'))HMt")) 

- (HMt"mix{t')))] , 
where the average value of an observable A is given by 



(9) 



{A) = 



N 



2^]^ cosh I 



fc=i 





Tr 




Wbt) 







After tracing the reservoir degrees of freedom from Eq. JHI, the influence functional acquires the form, 



T\x, y] = l-j^J^ dt' "^^''Yyk y) COS [ojkit' - t")] - ig{x, y) tanh [l^j sin [, 



UJk{t'~t")] 



(10) 



where 



/(x, y) = x{t')x{t") + y{t")y{t') - y{t')x{t") - y{t")x{t'), 
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and 

g(x, y) = x{t')x{t") - y{t")y{t') + y{t')x{t") - y{t")x{t'). 

Equation may be simplified by introducing a set of coordinates corresponding to the particle center of mass q 
and relative coordinate ^. Therefore, x — q + ^/2, y = q — ^/2, and the infiuence functional reads 

nx,y] = ^*"E"^' [m^in cosK(i' -t")] -2zg(<')e(i")tanh (^^^ sin - • (11) 

In the majority of cases considered we do not have enough information about the system which would allow us to 
perform the above summation in k. In order to overcome this difficulty, one usually introduces a phenomenological 
spectral density function, which correctly describes the bath in some well-known limiti^ Because we intend to follow 
this strategy, we will first calculate the form of the spectral density function associated to the two-level system reservoir 
and will replace it by some phenomenological guess afterwards. 

The spectral function can be obtained from 

J{uj) = Im id {-te{t - t') [F{t), F{t')])) , (12) 

where F{t) is the force produced by the particle over the thermal bath , 9(i) is the usual step function, and ^ stands 
for the Fourier transform. From Eq. I^J it is straightforward to demonstrate that F{t) = Ylk'^k^xk- Substituting 
then the expression for the force into Eq. II12II and using Eq. ijSl we obtain 

N / h \ 



fc=i 



Notice that the above expression is completely different from the oscillator spectral function. Here, all modes with 
energy well below ksT contribute less effectively to the spectral density than the "unoccupied" modes above the 
thermal energy. Therefore, any attempt to replace ifT^jl by a phenomenological guess should preserve this property. 

If we now replace Eq. lH^^jl into Eq. Illl|l , we obtain the influence functional 




I^^t^lll (coth ( i{t' W) COS [uj{t' - t")] - 2iq{t')i{t") sin \uj{t' - t")] 
TT 1^ \2kBT J 



where we have re-exponentiated the second order expan- 
sion used in Jsj. 

Next, we substitute Eq. (fTljl into Eq. 10 to obtain the 
superpropagator for the particle of interest, namely 



J 



I Vq exp<i -S'e//[g,e] } x 



exp < -- 



Jo 



$(t' - t")^{t')^{t")dt'dt" 



In the above equation, the effective action is given by 



and 



/ h \ 

J{lu, T) cos[w(t' - t")] coth I ) duj. 

\2kBT J 



After having traced the environment coordinates, we 
can derive from the effective action ltT5|l an equation of 
motion for the time evolution of the particle center of 
mass q and for the width ^ of the wave packet associated 
to it, 

Mq{t)+ [ K{t-t')q{t')dt' = eE{t), (16) 



where 



dt' 



Mq{t')i{i!)+ei{t')E{t') (15) 



+ / K,{t' ^t")q{t')at")dt"], 
Jo 



Ai-- 



J(w,T) sm[u; {t' -t")]du;, 



M^{t) + A{t- t')i(t')dt' = 



(17) 



where we have performed an integration by parts in order 
to explicitly show the viscous force of the TLSs reservoir 
acting on the particle and 

2 1"^ 

K{t-t') = - uj-'^J{uj,T)cos[Lj{t-t')]du;. (18) 
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The last step consists in solving Eqs. Hl(i|l and II17|I . 
However, we must first specify the form of the spectral 
density. Because we intend to keep the problem as gen- 
eral as possible, we will choose a form for the spectral 
density which retains the functional T-dependence given 
by Eq. ((T3|l . We then assume that 



(19) 



where 17 is the cutoff frequency already introduced in the 
discussion of the model, 7 is a constant defining the par- 
ticle coupling strength to the TLSs and s is a c-number. 
The reader must be warned that 7 used in this paper is 
actually Aj/ir where 7 is the usual relaxation constant of 
the particle motion. 

Expression (jl9ll retains the main properties of the func- 
tional form Ijl3p . namely the fact that the temperature 
determines which are the statistically relevant modes. 
Therefore, for a given value of and depending on the 
temperature, the particle may simply have no interaction, 
on average, with the reservoir. Notice that this choice is 
at variance with the one employed in Ref. whre the 
author maps the system onto a bath of oscillators. 

By replacing Eqs. Ijl8|l and (jl9ll into lll(i|l . we obtain the 
following equation for the velocity of the particle center 
of mass: 



v{t)+^ / r{t-t')vit')dt' = 



M ' 



(20) 



with the damping function 



r{t-t') 



tanh 



huj 
2kBT 



cos[uj{t — t')] du. 



(21) 

As it can be observed from Eq. II2()|I . after tracing the 
two-level system reservoir coordinates, we have obtained 
an equation of motion for the particle center of mass 
in which the thermal bath has the same effect as that 
of a viscous fluid. It should be noticed that for s = 1 
and zero temperature, Eq. II21II reduces to the oscillator- 
bath damping function, which is memoryless in the limit 

— > 00 or, in other words, in the long time regime. 
Indeed, in this case 



r(t-t') 



lim 



C08{uj[t - t'])dLO = TT6{t - <'), 



indicating that the damping is a purely instantaneous 
function. 

However, if s 1 we realize that even at zero tem- 
perature and with the cutoff frequency going to infinity, 
it is impossible to obtain a damping function without 
memory. In this situation the problem becomes non- 
Markovian and the particle dynamics can not be de- 
scribed in terms of damping and diffusion coefficients. 
In the following, we will investigate the transport prop- 
erties of a system of non-interacting particles described 
by Eqs. lf2?Hl and 1I2TI1 through the evaluation of its optical 
conductivity. 



IV. THE OPTICAL CONDUCTIVITY 

The current associated to distinguishable non- 
interacting particles described by the equation of motion 
l|2(11l reads j = ev{t), where j satisfy 



(22) 



This equation may be solved by the well-known method 
of the Laplace transform. If we assume that initially 
there is no current in the system, Eq. Ij22|l may be written 

as j{z) = a{z)E{z), with 



aiz) 



M[z + jr{z)Y 



where r(z) is the Laplace transform of T{t — t'). For all 
allowed values of s, T{z) can be written as (see appendix) 



T{z) = : \; J^' tan ' 



(s + l)z2 



\2kT 



(23) 



where F denotes hyper-geometric functions given by 

F{s,x) = 2Fl{l,^^,^^,x) and A„ = {2n ~ lykeT /h, 
with n e N. 

The optical conductivity may now be promptly ob- 
tained by substituting lf23jl into 



alu}) = lim Re 



eVM 




z + ^T{z) 


z—5—iuj 



(24) 



Although Eqs. Il23ll - lj24|l allow us to compute <T(aj) for 
any s > 0, it may be enhghtening to consider the prob- 
lem for specific values of s, for which the hyper-geometric 
series in II23II converge to simple functions. In this way we 
will be able to proceed analytically in the investigation 
of the optical conductivity and obtain, if there are such 
contributions, the temperature dependence of the Drude 
weight and the incoherent conductivity in the whole fre- 
quency range. In the following we will focus on particular 
cases illustrating the super-ohmic (s > 1), ohmic (s = 1) 
and sub-ohmic (s < 1) situations. 



1. Super-ohmic case, s = 2 

In this specific case, the hyper-geometric functions in- 
volved in Eq. (I23II acquire a simple form 



2Fi(l,3/2,5/2,-a;^) = - -- 
X \ x 



1 arctan x 
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Figure 2: (a) Drude weight as a function of T for (1 = 1, the 
continuous, dotted and dashed lines correspond to 7 = 0.2, 
7 = 0.8 and 7 = 1.6 respectively, (b) Temperature behav- 
ior of — l/lncr^'^, the continuous, dotted and dashed lines 
correspond to 7 = 0.2, 7 = 0.4 and 7 = 0.6 respectively, 
(c) — InT for low temperatures and 7 — 0.2. In 

all cases it is assumed ft = fc_g = 1 and the Drude weight is 
measured in units of ctq. 



and the Laplace transform of the damping function re- 
duces to 



r(.) 



4^ °° (2n — 1) arctan 



no, 



irfcBT(2ri-l) 



n=l 



(2--l)^-(^)' 



,T2\ / hz 
-z arctan — tan 



(25) 



Substituting then Eq. ll25ll into lf2ijl . we obtain the fol- 
lowing form for the optical conductivity: 



(26) 



where {T) represents the Drude weight, whereas 
a''"'^{uj,T) stands for the incoherent contribution. Let 
us first analyze (T), which is given by 



O (J j = CTo 



1^ °° arctan 

TT ^ 



{2n- 1) 



(27) 

where ao — 7re^/M is the Drude weight of the free parti- 
cle. Although this sum can not be performed exactly for 
all temperatures, we may derive an analytical expression 
for it in the regime ksT ^ Hfl. We then find 



nw 



(T) 



4jhn \ " 1 

ir'^kBT/^ (2n-l)^ 



(l+75l^) 



-, (28) 



where the damping term provides a small correction to 
the particle mass. Notice that Eq. lf28jl correctly repro- 
duces the free particle behavior at high temperatures. 
This result is in agreement with the fact that as the tem- 
perature is raised the high frequency TLSs will play a 
major role in the composition of the spectral function 
which justifies the mass correction as due to the adia- 
batic approximation. In other words, the high frequency 
TLSs adiabatically dress the moving particle. 

A numerical evaluation of the sum involved in Eq. lf27jl 
yields the general behavior of the Drude weight as a func- 
tion of temperature (see Fig. El^a)). Observe that the 
conductivity grows with temperature until it saturates 
at the value of a free particle in the T cx) limit. In 
order to investigate the functional growth of the conduc- 
tivity with temperature, we have plotted —l/\na^^{T) 
vs T in Fig. Efb). The plot is linear in T, except in 
the low temperature region. We may therefore conclude 
that for high temperatures the Drude weight behaves as 
oc exp[— l/(a -I- (3T)], where the constants a and 
/3 can be determined from the plots and depend on the 
values of and 7. On the other hand, the vanishing of 
{T) as the temperature is lowered was already ex- 
pected from Ij27|l . because for strictly zero temperature 
the sum to be performed is positive and divergent. 

The functional reduction of the conductivity as the 
temperature decreases may be determined by consider- 
ing a certain N — 2ri ~ 1 such that hfl/irkBTN ^ 1. In 
the limit of T ^ the value of N will be of the order 
of hfl/kBT. In this situation Eq. lf27jl may be approxi- 
mately written as 



„DW /rr.^ 



a-tt; °° arctan 



7rfci3T(2n-l) 



2n- 1) 



where 



5 = E^="+2^'"H^+V2), 



(29) 



(30) 



n=l 



with a = (C + ln4)/2, 



m ^ 

C= lim V- 

TT). — ^DO ' ^ h 



In TO, 



fe=i 



and ^^^^ denote the polygamma function of zeroth order 
given by ■(/'("^(a;) = dxlnr{x). In the zero temperature 
limit the sum given by Eq. H3()|l dominates the behavior of 
Eq. lf29jl . Using then the Stirling expansion and assuming 
that n is of the order of Ml/ksT, we finally obtain that 
the Drude weight in the low temperature limit behaves 
as 



nw 



(T^O) 



M 



(31) 



In the expression above £ is a numerical factor that can 
determined from the plot of l/cr^'^vs — l/lnT for low 
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temperatures. Inspection of Fig. I^Jc) indicates that the 
value of £ is approximately 1. Therefore we can conclude 
that the Drude weight for low temperatures behaves as 



nw 



(r->o) 



( hO. \ 



(32) 



In 



Actually the effect of finite 7 can even wash the Drude 
weight out as in the nai've Drude model for electric con- 
ductivity of metals. Nevertheless, in our specific model, 
although is reduced, it is still finite even in the 

presence of damping for finite T. Moreover, as the ra- 
tio hfl/kBT is directly proportional to the number of 
TLSs in the lowest energy state, it is expected that as 
the temperature is lowered this number rises, increasing 
the ability of the particle to lose energy and consequently 
leading to a smaller value of . 

Now we turn our attention to the behavior of the in- 
coherent part of the optical conductivity, which is given 

by 



a [y>) = ctq- 



7 tanh 



2/cbT 



ein-Lo) 



^ tanh (2!^^ 



with 



=1 



(2n — 1) arctan 



7rfcBT(2n-l) 



TT 



7 

tanh 



(2n 
fku 



1)2 + {hujI-KkBTf 



(33) 



(34) 



In 



The presence of the step function in Eq. Ij33|l ensures 
that (T*"'^(a;) is exactly zero above the cutoff frequency. 
Although for frequencies of the order of SI our approach 
may provide non accurate results, the zero conductivity 
in the > SI region agrees with the fact that all particle 
transitions between states with energy difference larger 
than Ml are forbidden. This effect is a consequence of 
having limited the thermal bath phase space by an abrupt 
cutoff frequency. 

In order to discuss the main features of the incoherent 
part of the optical conductivity given by Eqs. ll33|l - l(3ijl 
we begin by calculating the temperature dependence of 
the dc conductivity adc = <j^^^(uj = 0), namely 



crdc{T) _ h'j 



CTO 



47=57 °° arctan 

1 H -Y^ -, — 

TT ^ (2n 

n=l ^ 



TTkBT(2n-l) 



1 



The high temperature Hmit of this expression can be writ- 
ten as 



CTO 



1 



?t~0 



2-KkBT 




Figure 3: (a) ct'"'^ as a function of u) for different tempera- 
tures with 7 — 0.2 and (1 = 1. The continuous, dashed and 
dotted lines correspond to T = 0.09, T = 0.5 and T — 2.0 
respectively, (b) Details of cr(a;) near uj = ft. (c) o"*"'^ as a 
function of to for different coupling strengths. The continu- 
ous, dashed and dotted lines correspond to 7 = 0.4, 7 — 1.5 
and 7 = 3.5 respectively, with T = 0.02 and Q = 1. (d) a'"" 
vs a; for 7 = 0.2 and T = 0.05 for different cutoff frequen- 
cies, the continuous, dashed and dotted lines corresponds to 
r2 = 1, $1 = 1.5 and = 2 respectively. In all cases the opti- 
cal conductivity is measured in units of ctq and it is assumed 
h — ks = 1. 



which correctly goes to zero as the temperature incre- 
ses in agreement with the free particle behavior. The 
interesting point comes from the divergence at zero tem- 
perature. This singular behavior has been observed in 
classical non-integrable nonlinear systems where the cur- 
rent correlation decays to zero in the long time limit, but 
so slowly, that the integral over time, which yields the dc 
conductivity, divergesii^ In the present case this result 
can be understood by reaHzing that the long time behav- 
ior of the damping function II21|I involves low frequency 
modes, but the latter do not have enough spectral weight 
for the s — 2 case to make the current decay. 

A further understanding of the problem can be 
achieved after the analysis of Fig. where the fre- 
quency dependence of tT™^(a;) is plotted for some par- 
ticular cases. Notice in Fig.EJa) that as the temperature 
increases, the conductivity is reduced in the entire fre- 
quency range. This is in agreement with the fact that 
as one moves toward the limit in which fcsT Ml the 
particle approaches the free behavior and eventually the 
incoherent conductivity vanishes, see Eqs. Ij33|l and II34II . 
retrieving the free particle result. 

Let us now inspect the behavior of (t™^([x)) when uj 
n. Fig. Etb) shows a zoom of Fig. Ufa) for lu around 
ri. Near the cutoff frequency the conductivity reaches 
a maximum and then continuously falls to zero as one 
approaches fl from below. This non-monotonic behav- 
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ior of the optical conductivity as a function of the fre- 
quency can be basically attributed to memory effects of 
the damping function lf2T|l . This effect is responsible for 
the dephasing of the different contributions of the bath 
reaction over the particle. As a matter of fact, if we rean- 
alyze Eq. l(33ll . we observe that the classical impedance 
Z{uj,T) of the reservoir (for uj < il) corresponds to 



Z = ag-i 



2cjQ2(a>,r) 



7 tanh 



\2kBT J 



— - — tanh f 



Although this expression can not be clearly separated 
into reactive and resistive contributions, one of the com- 
ponents of the former will be 



<Y(w,r) oc 



2t^ Q2(w,T) 



7 tanh 



where Q(a;,T) is given by Ij^MII . From the definition of 
Q(a;, T), we observe that part of the reaction of the bath 
over the particle is nothing but a competition of terms 
with a TT phase difference between them, which causes 
this kind of resonance-like behavior. The effect described 
above does not appear in the oscillator model, where the 
bath reacts as a pure resistance due to its memoryless 
character, and leads to a monotonic behavior of a^'^'^iu). 
It is also expected that as the coupling strength of the 
particle- reservoir interaction increases, the conductivity 
will decrease. This effect is illustrated in Fig. |3Ic), where 
it can be seen also that for strong enough coupHng the 
optical conductivity becomes a monotonic function. It 
remains to mention that for a fixed value of the tempera- 
ture, an increase in the value of the frequency cutoff leads 
to a displacement of the conductivity edge to higher fre- 
quencies and to a decrease of cr™'^ in the whole frequency 
range. This effect is illustrated in Fig. El^d) and it is a 
consequence of the higher number of states to scatter the 
particle, as the value of the cutoff frequency is increased 
at a given temperature. 

All the previous analysis developed for the s = 2 spe- 
cific case illustrates to some extent the transport prop- 
erties of the system in the super-ohmic regime (s > 1). 
However, it is worth mentioning that there are signifi- 
cant differences in the low frequency region of the optical 
conductivity as one moves from l<s<2tos>2 sit- 
uation - see Fig. 31 for instance. The divergence in the 
dc conductivity is suppressed for s > 2 while it persists 
for s < 2. This is again a consequence of the different 
spectral weights for the low frequency modes as the value 
of s is changed. 

To conclude this part of the analysis we summarize our 
main findings. In the specific case of s = 2 the system 
behaves as a perfect conductor at T = because of the 
infinite Odc conductivity. However, this fact is not re- 
fiected in the value of the Drude weight, which goes to 
zero as T decreases. This contradiction in the classifi- 
cation into metallic or insulator states according to the 




Figure 4: cr'"'^ as a function of - in units of ctq - for differ- 
ent values of s in the super-ohmic regime. Notice that as s 
increases, the value of the Odc decreases for a given frequency. 
All graphics were generated assuming 7 = 0.2, T = 0.09, 
n = 1 and h = kB = 1. 



value of the Drude weight is only apparent. This is so 
because the current carriers in our approach are not in an 
eigenstate of the system invalidating the characterization 
of conductors and insulators in terms of the behavior of 
fjDW Another point to be mentioned is the fact that for 
s 7^ 2 the dc conductivity, at T = 0, strongly depends on 
the s value, see Fig. 31 Finally, above the cutoff fl, the 
optical conductivity is always zero at any temperature, 
independent of the s value. 



2. Ohmic case, s = 1 

This case is of particular importance because it allows 
us to explicitly show the differences with the oscillator 
model discussed in Refs. j3| and In this specific 
situation it is useful to notice that 



In a; 



2Fi(l,l,2,-x2 



Therefore, the general expression ^'21M for the Laplace 
transform of the damping function becomes 



r(z) 



1 / 








- tan 




H 




2 ^ 


\2kBT) 







(35) 



2hz 
n'^ksT 



°o In 1 1 



[hn/{2n-l)TTkBTf^ 



(2n- 1)2 - {hz/TxkBTf 



By substituting l(35|l into 1I2H1 . we obtain the opti- 
cal conductivity as a sum of a coherent part given by 
a^^{T)5{uj), where 



-rDW 



ln<^ 1 



{2n-l)iTkBT 



\ 



(2n- 1)2 



(36) 
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and an incoherent contribution of the form 

(77r/2) tanh (5!^) (^o^{^ - uj) 



(a) 



where 

n = uj- 



7^ tanh (2^) /2]2+7l2(u;,T) 



(37) 



^ tanh 



2 



E 



hi|l - 



hi 



(38) 



(fif7/(2n- l)7rfcBT)^| 



(2n- 1)2 + {hiu/TTkBTy 



Once again the Drude weight can not be obtained ex- 
actly for all temperatures, except for ksT 3> hfl. In this 
limit we find 



nw 



(T) = ao 



2 — 



96rJ KkeT 



(39) 



which correctly reproduces the expected free particle con- 
ductivity in the infinite temperature regime, where the 
incoherent contribution goes to zero [see Eqs. (jH7|l and 
Il88p ]. If one compares Eqs. Ij28|l and Ij39|l . one realizes 
that for the same temperature the Drude weight in the 
s = 2 case is lower than in the s = 1 . This is nothing but 
the effect of the stronger coupling between the particle 
and the reservoir modes as value of s -characterizing the 
thermal bath- is increased. 

The low temperature properties of the Drude weight 
can be derived from Ij36|l by considering a value of n = n 
such that hfl/kBT 3> n. In this case the sum involved in 
II36II can be written approximately as 



S 



1 



ln7r(2n - 1) 



V /csT / (2n- 1)2 ^ (2n-l)2 

^ 71—1 ^ ' n— 1 ^ ^ 



(2n- 1)2 



Notice that for T ^ the value of n goes to infinity and 
therefore the first term in the rhs of Ij4fl|l dominates the 
sum. Substituting this result into l(3fi|l we obtain that for 
low temperatures the Drude weight goes to zero as 



nw 



(T 0) oc 



keT 



hrf In 



no 



fesT 



(41) 



This expression correctly reproduces the well known re- 
sult for the oscillator model;^ namely there is no coherent 
contribution to the conductivity when goes to infinity. 
This result could have been foreseen from the analysis of 
the particle dynamics [see Eq. Ij21ll for s — \ and T — > 0] . 
It should be stressed that for s — \ goes to zero 

faster than for s = 2 [see Eq. l(32)l ]. as a consequence 
of the stronger interaction of the particle with the low 
energy modes. 




5 
T 



10 





(d) 



I -. I 1 

1 

— '.X 


1 









0.5 

0) 



Figure 5: (a) rj^^ as a function of T for different coupling 
strength with 0,-1. The continuous, dotted and dashed 
lines correspond to 7 = 0.2, 7 = 2 and 7 = 5 respectively, 
(b) a{Lo,T) vs uj for different temperatures with 7 = 0.2 and 
n = 1. The continuous, dotted and dashed lines correspond to 
T = 0.001, T = 0.02 and T = 0.09 respectively, (c) The same 
as (b) assuming Q. = 100. (d) a"''' vs ctJ for = 1 and T = 
0.02 for different coupling strength, the continuous, dashed 
and dotted lines corresponds to 7 = 1, 7 = 0.2 and 7 = 0.09 
respectively. In all cases the conductivity is measured in units 
of (70 and it is assumed /i = fcs = 1. 



After having examined the temperature dependence of 
the Drude weight for ohmic dissipation in the analyti- 
cally accessible limits, we proceed to a numerical eval- 
uation of Eq. l(3fijl . Fig. Ela) shows the behavior of 
the Drude weight as a function of temperature for differ- 
ent coupling strength values. Observe that for any given 
temperature a^"^ decreases as the interaction strength 
becomes stronger. Indeed, for stronger interactions the 
momentum transferred to the reservoir is larger and con- 
sequently the conductivity is reduced. This also implies 
that the free particle behavior will be reached at higher 
temperatures as 7 increases. 

Our next step is the analysis of the incoherent part 
of the optical conductivity. The first point to be men- 
tioned is the infinite dc conductivity of the system at 
any finite temperature. This result is obtained taking 
the limit of zero frequency in the expressions Ij37|l and 
ll38)l and can be explained by the fact that the dc con- 
ductivity is determined by the low frequency modes in 
the damping function II21|I . We have plotted in Fig. EJb) 
cr*"^ vs uj for different temperatures with a finite cut- 
off frequency = 1. The forbidden particle transitions 
involving energy exchange above hSl leads to zero con- 
ductivity for > f2, as already obtained in the s — 2 
case. One may also notice that as the temperature in- 
creases the conductivity grows for small values of uj. This 
can be understood by recalling that as the temperature 
rises, the particle-reservoir interaction - which can ac- 
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tually change the particle momentum - is less eflective 
for the low energy modes and consequently the conduc- 
tivity also rises for such frequencies. It should be also 
mentioned that as the temperature is lowered, cr™^((jj) 
acquires a Lorentzian form. This behavior becomes ev- 
ident in Fig. Etc), in which we considered = 100. At 
lower temperatures one approaches the oscillator model. 
Indeed, taking the limits ^l ^ co and T — > 0, we obtain 



7r7/2 



(7^7/2)- 



■ CTo 



which is the well known result for the oscillator model. 
Finally, Fig. E^d) shows the effect of varying the coupling 
strength. The main feature to be observed is the reduc- 
tion of the conductivity as 7 increases, consistent with 
the fact that it is then easier for the particle to transfer 
momentum to the reservoir. 



3. Sub-ohmic case, s = 

Another interesting result, which differs considerably 
from the cases already discussed, arises when s = 0. In 
this case, the hyper-geometric function reads 



2Fi(l,l/2,3/2,-x2) = 



arctan x 



and the Laplace transform of the damping function given 
by Eq. (l23|l acquires the form 



A' — z~ 



rctan (^) 



arctan 



(ft) 



An 



(42) 

(2n - l)kBTn/h with n 



where we have used that A„ 
integer. 

Substituting now Eq. l(l2|l into lf2ijl . we obtain an 
expression for the optical conductivity which has zero 
Drude weight at all temperatures and an incoherent part 
of the form 



inc / \ _ 



with 



27u;tanh (arfr) ^(^ " ^) 
g^{uj,T) + [7r7tanh {huj/2kBT)]' 



8h^jLu^j2 arctan (nf^/iVTrfesT) 



(43) 



7r3fc2r2 N [N^^ + [hw/TTkBTf 



-I-7 tanh 



fku 
2kBT 



In- 



\VL-uj\ 



\VL- 



(44) 



The dc conductivity, obtained by taking the ^ 
limit in Eqs. Ij43ll - lj44|l . may be written as Cdc/fo — 
AkBT /ir'^tvy. From this result we conclude that the sys- 
tem behaves as an insulator with zero dc conductivity at 
T = 0. 




Figure 6: (a) Incoherent contribution to the optical conduc- 
tivity for different temperatures with 7 = 0.2 and ^ — 1. The 
continuous, dotted and dashed lines correspond to T = 0.05, 
T = 0.2 and T = 0.4 respectively, (b) Detail of (a) near 
the cutoff frequency, (c) cr'"'^ vs co for different values of 
il. The continuous, dashed and dotted lines corresponds to 
fl = 100, = 1.5 and = 1 respectively, in all cases 7 = 0.2 
and T — 0.05. (d) cr*"'^ vs lj for different values of the cou- 
pling strength. The continuous, dashed and dotted lines cor- 
responds to 7 = 0.2, 7 = 0.25 and 7 = 0.3 respectively, with 
n = 1 and T = 0.09. In all cases the optical conductivity is 
measured in units of o"o and it is assumed h = ks = 1- 



Inspection of Eq. l(i3jl yields that the incoherent part of 
the conductivity is finite only below the cutoff frequency, 
just as in the cases discussed previously. Another impor- 
tant feature is the way in which one may recover the free 
particle behavior. In order to analyze this limit it should 
be noticed that at high temperatures Eqs. H43|l and Ij44|l 
may be written as 



= (^) 



c 



M w2 + (^2 ' 



(45) 



where C = irhj/AkBT. In the limit C ^ 0, which corre- 
sponds to infinite temperature or zero coupling strength, 
Eq. Ij45|l goes to {Tre^ /M)5(uj), correctly reproducing the 
optical conductivity of a free particle. 

The general behavior of (t™'^(w) is displayed in Fig.EI 
for different values of temperature, coupling strength and 
cutoff frequency. Fig. Eta) shows that as the tempera- 
ture increases, one approaches the frequency dependence 
given by Eq. I^IJ, which has a maximum at w = 0. In this 
regime, the classical impedance of the reservoir is purely 
resistive, analogous to the oscillator-bath model. How- 
ever, as the temperature decreases, the value of Z{uj) is 
determined by a competition between terms completely 
out of phase, see Eq. II44|I . reaching the maximum value 
for u) ^ 0. For T ^ 0, the frequency at which the 
conductivity exhibits a maximum becomes closer to the 
cutoff frequency, characterizing an extreme non-Drude 
behavior. As one approaches fl from below, no matter 
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the value of the temperature, the conductivity smoothly 
falls to zero. This effect, which is related to the lack 
of available states in the reservoir to scatter the parti- 
cle, is shown in Fig. E^b). It is also worth mentioning 
the small variations of the maximum conductivity as the 
cutoff frequency is changed, as illustrated in Fig. EIc). 
Finally, Fig. Eld) shows the effect of varying the value 
of the coupling strength 7 for a given temperature and 
cutoff frequency. An increase in the interaction strength 
leads to a reduction of the conductivity. This is an ex- 
pected result because for larger values of 7, the particle 
can transfer momentum to the reservoir in an effective 
way, and as a consequence, the conductivity is reduced. 

At this point we can conclude the s = analysis point- 
ing out that in this situation the conductivity is always 
incoherent and non zero just below the cut off frequency 
with a dc value that vanishes as T goes to zero. In order 
to verify if those features are common to the entire sub- 
ohmic regime, we have numerically evaluated the com- 
bined Eqs. Ij42|l and Ij24|l for different values of s < 1, see 
Fig. CJa). As it can be observed, the general behavior 
of the conductivity resembles the s = case (see Fig. 
6), except in the low frequency region. This fundamental 
difference is illustrated in Fig. Clb) showing that cr™'^ is 
singular at oj — 0, although our numerical calculation in- 
dicates that the Drude weight remains zero at any finite 
temperature and value of s. In order to investigate in 
more detail the general behavior of ct™'^ in the low fre- 
quency region, it is useful to notice that near oj = the 
optical conductivity does not depend on the value of the 
cut off frequency, see Fig.CJc). Therefore, we may simply 
take the limit ^ 00 in the Laplace transform of the 
damping function given by expression ll23|l and obtain an 
expression for T{z) for any < s < 1, namely 



r(z) 



2TrkBT z 



hC0s{TTs/2) 



00 



A? 



(46) 



Substituting the expression above in Ij24|l we obtain 
the optical conductivity, which has a zero Drude weight, 
besides an incoherent contribution given by 
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Figure 7: (a) Incoherent contribution to the optical conduc- 
tivity in the sub-ohmic regime, in all cases it was assumed 
7 = 0.2, T = 0.09 and = 1. (b) Detailed form of the con- 
ductivity in (a) near uj = for s = 0.6 and s = 0.3. (c) u'"'^ 
vs UJ for different values of fi, in all cases T — 0.09, 7 — 0.2 
and s = 0.3. The optical conductivity is always measured in 
units of (70 and it is assumed h = ks = 1- 



= (w - 0) _ 4fcsT cos2(7rs/2) 



It is worth mentioning that this expression is also valid 
for the s = case (see Eq. l(l2|l . for instance), correctly 
reproducing the finite adc{T) value already found. 

With the analysis of the — > 00 limit we conclude that 
the finite dc conductivity obtained for the s = case is 
an exception within the sub-ohmic regime, in which the 
dc conductivity is, in general, divergent. This discussion 
summarizes our main findings in the study of the optical 
conductivity for cases in which < s < 1. 



CONCLUSIONS 



= (u;) 



27^^-3 tanh(5^) 



AlC 



TTjuj'^ 1 tanh( 



2k bT 



where 



IC = 1 



— tan( — )tj'*"nanh — 

2 ^ 2 ' \2kBT 



E 



{2n-lY 



2--f{TTkBTf 

'^^cos(f) ^^i2n-lY + {^Y 



(47) 



(48) 



This result is in complete agreement with the previ- 
ous numerical analysis for finite cutoff frequency and 
promptly allows us to derive the way in which the con- 
ductivity diverges as one approaches w = 0, 



We have studied the transport properties of non- 
interacting particles coupled to a set of two-level systems 
described by a temperature-dependent spectral func- 
tion. Our approach was based on the Feynman- Vernon 
functional-integral formaHsm which allows us to trace out 
the two-level system coordinates and obtain an effective 
equation of motion for the particle in terms of the phe- 
nomenological TLSs spectral function. 

The evaluation of the transport properties were per- 
formed considering different values of the power s of 
the generic spectral function, which is proportional to 
ijj^ . Two different regimes were found; the super-ohmic 
(s > 1) and ohmic (s = 1) in which the optical conductiv- 
ity has both, coherent and incoherent contributions and 
the sub-ohmic (s < 1) where the conductivity is strictly 
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DW 


DW 


inc 


inc 




Ode 






Ot~>0 




Super-Ohmic 


diverges for T — > 0" 


zero 




Anomalous insulator - non Drude 


zero 


Ohmic 


diverges for finite 


zero 


cro 


Insulator - Lorentzian form'' 


zero 


Sub-Ohmic 


diverges for finite 






Perfect insulator - non Drude 


a-o5(w) 



"For s < 2, in other cases it is finite 

''For Q ^ oo reproduces the oscillator results 

■^Except in the s=0 case, where the o-^c goes to zero 



incoherent. In the first case the Drude weight goes to 
zero as the temperature decreases in a logarithmic fash- 
ion and the incoherent part shows a strong non-Drude 
behavior with a divergent dc conductivity for T = 0, ex- 
cept for the ohmic case in which the zero temperature 
dc conductivity is finite. In the entire sub-ohmic regime, 
the Drude weight is always zero at all temperatures. In 
contrast, the dc conductivity is only finite in the s = 
specific case, going linearly to zero as the temperature 
decreases, whereas in general it is divergent at all tem- 
peratures. It should be mentioned that in this regime, 
the general behavior of cr(a;) is strongly non-Drude, with 
the highest conductivity arising at finite values of lu. All 
these properties strongly differ from the simple behav- 
ior induced by an oscillator thermal bath on a system 
of non-interacting particles. The latter exhibits no tem- 
perature dependence and zero Drude weight with a finite 



dc conductivity. All this results are summarized in the 
table above. Concerning the limitation of the approach 
employed here, it is important to stress that it is accu- 
rate for long times (or low frequencies) compared with 
the response time of the system, which is basically . 
Therefore in all cases where the cutoff frequency was as- 
sumed to be finite, the values of the optical conductivity 
near oj — ^ should be interpreted carefully. 

Finally, it is worth pointing out that the model an- 
alyzed here could be relevant for the study of the dy- 
namics of a particle in the presence of a distribution of 
locally quartic plus quadratic potentials. We expect that 
this simple model may contribute to the understanding 
of the transport properties of low dimensional systems at 
low temperatures, where the physics is often dominated 
by defects and impurities. 
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Appendix 

In order to calculate the Laplace transform of T{t) we 
will use the expansioni^ 



I Tix\ Ax s—^ 1 



2 7 TT ^ (271 - 

n=l ^ 

which allows us to rewrite T{z) as 



Viz) 



^I(n,f]), 



(49) 



(^'2+^2)(^,2 + ^2)^ 



with A„ = (2n— l)nkBT/h and n e N. This integral can 
be split into two terms. 



1 



The new integrals are easily performed, yielding 



E 



(-ir 



^ — ' 2m + s — 1 \ u 

m=o ^ 



We now rewrite the expression for I(n, Q,) as 



^ oc 



^ oo 
72 S 



(-ir 



z'^ z — ^ 2m + 5 — 1 \ z 

m— 1 



n 



i-iy 



n 



" m=l ^ 



2m-2 



(50) 



n=0 



Substituting this equation into ijifljl and expressing the 
sums involved in Ij5()|l in terms of the hyper-geometric 
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functions we obtain 



T{z) 



+ 1) t{ 
2i^l(l, 



2 ' 2 ' 



By performing the sum in the first term of Eq. HR1|I and 
introducing the notation 2^i(l, "4^, ^3^, a:) — F{s,x), 
we finally obtain Eq. lf23jl . 



1 + 6- 3+6- 

2 ' 2 ' 



\ 2 
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